EFI-95-65 
FERMILAB-PUB-95/345-T 
ANL-HEP-PR-95-80 
WIS-95/50-Oct.-Ph 
hep-ph/9511363 
November 1995 



ANGULAR DISTRIBUTIONS AND LIFETIME DIFFERENCES 

^ : IN B,^ J/^lj<i) DECAYS [| 

a^ : 

Amol S. Dighe f\ 
Q I Isard Dunietz q| 



> 

cn 
cn 



Harry J. Lipkin^ f\ 
and 

Jonathan L. Rosner ^' ^ 
ABSTRACT 



The strange B meson Bg = bs and its charge-conjugate Bg = bs are ex- 

^ I pected to mix with one another in such a way that the mass eigenstates B^ 

■ ( "heavy" ) and B^ ( "hght" ) may have a perceptible hfetime difference of up 

p ■ to 40%, with the CP-even eigenstate being shorter-hved. A simple transver- 

p . sity analysis permits one to separate the CP-even and CP-odd components 

D ! of -Bs — > J/ipcf)-, and thus to determine the lifetime difference. The utility of 

^ ' a similar analysis for B^ J/ipK*'^ is noted. 



X 

^ I The Cabibbo-Kobayashi-Maskawa picture of weak charge-changing transitions |1| 

predicts the strange B meson Bg = bs and its charge-conjugate Bg = bs to mix with one 
another with a large amplitude. The mass eig enstates 5f ("heavy") and B^ ("light") 
with masses m{B^) = mn and m{Bg) = nii are expected to be split by Am = ruH — 
ruL ~ 25f , give or take a factor of two |[, where f = (F^ + Tl)/2 ^ T{B°) (5° = bd) 
and Th,l = ^{B^ , Bg). The measurement of such a large mass difference poses an 
experimental challenge. 

To a good approximation, CP violation can be neglected in calculating the mass 
eigenstates, in which case they correspond to those -B^^-* of even and odd CP, with 
Bg = -B]'*'^ and B^ = B^~^ as we shall see. The decay of a Bg meson via the quark 



-'^To be submitted to Physics Letters B 

^Enrico Fermi Institute and Department of Physics, University of Chicago, Chicago, IL 60637 
•^Theoretical Physics Division, Fermi National Accelerator Laboratory, Batavia, IL 60510 
"'High Energy Physics Division, Argonne National Laboratory, Argonne, IL 60439 
^Department of Particle Physics, Weizmann Institute of Science, Rehovoth, Israel 



1 



subprocess b{s) ccs(s) gives rise to predominantly CP-even final states 0. Thus the 
CP-even eigenstate should have the greater decay rate. An explicit calculation M gives 



r(5W) - r(fi(-)) /• 



°-^^200 MeV)2 ' 

where fs^ is the Bg decay constant (in a normalization in which f^^ = 132 MeV). In 
one estimate fs^ = 225 ± 40 MeV, while a compilation of lattice results |^ obtains 
/b, = 201 ± 40 MeV (90% c.l. limits). The upper limit of 40% for is based on an 
estimate of the maximum possible contribution from the b{s) ccs(s) subprocess [|, 0]. 

The ratio of the mass splitting to the width difference of strange 5's is predicted to 
be large and independent of CKM matrix elements P, H] (to lowest order, neglecting 
QCD corrections which may be appreciable): 



Am _ 2 m^h{mj/M^) / 8 ml' 
AT ~ ~3^ V ~3^ 



-1 

I --200 , (2) 



where AF = Th — Pl. Here h{x) decreases monotonically from 1 at a; = to 1/4 as 
X — > oo; it is about 0.54 for mj = 180 GeV/c^. In view of the sign in Eq. (|1]) and 
since Am > by definition, we then identify = B^^^ and B^ = Bl~^ 0. If the 
mass difference Am turns out to be too large to measure at present because of the rapid 
frequency of Bg — Bg oscillations it entails, the width difference AP may be large enough 
to detect. The possibility of a value of AP/P for strange B mesons large enough to 
measure experimentally has been stressed previously 0, |^, |^ . 

One can measure f using semileptonic decays, while the decays to CP eigenstates 
can be measured by studying the correlations between the polarization states of the 
vector mesons in -B^^^ J / ip'P- (For similar methods applied to decays of other spinless 
mesons see, e.g., Ref. ||10[-) In the present note we describe a means by which the J/ijj4> 
final states of definite CP in Bg decays may be separated from one another using a sim- 
ple angular distribution based on a transversity variable |Tl|, |1^, |13|. This transversity 



variable allows one to directly separate the summed contribution of the even partial 
waves (S, D) from the odd one (P) by means of their opposite parities. The CDF Col- 
laboration has recently reported the first angular distribution analysis of the decay 
Bg J / ip(j), obtaining a separation into longitudinal and transverse helicity amplitudes 
without making a statement yet about the CP-even and CP-odd contributions. 

We summarize our main result. Consider the final state J /ipip — ^ i~^i~K^K^ , where 
i = e or fi. In the rest frame of the J/ip let the direction of the define the x axis. 
Let the plane of the K^K~ system define the y axis, with py{K^) > 0, so the normal 
to that plane defines the z axis. (We assume a right-handed coordinate system.) We 
define the angle 6 as the angle between the and the z axis. Then the time-dependent 
rate for the J/ip(f) mode is given by 

rf^F 3 _ /. n „\ 3 



■p(t) (l + cos^ 6*) + -m(t) sin^ 6 



dcosO dt 8 ^ ^4 
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= ^[p(t) + 2m(t)] + ^[p(t)-2m(t)]cos2^ , (3) 

where 

p(t) = p(0)e-^^* (CP even) , m(t) = m(0)e-^«* (CP odd) , (4) 

so that the probabihty of having a CP-even [CP-odd] state at proper time t is given by 
pit) /{pit) + m(t)) [m(t)/(p(t) + m(t))]. The angular distribution is normahzed in such 
a way that 

^^rf(cosg) , , =p(t)+m(t) . (5) 



(it 7-1 dcosO dt 

As t increases, one should see a growth of the sin^ 9 component. The angle 9 is an 
example of a transversity variable, whose utility for the determination of CP properties 



of multi-particle systems was pointed out some time ago ||T5 . 

The zero-angular-momentum states of two massive neutral vector mesons such as 
J/ip and 0, both with the same CP (in this case, even) consist of two with even CP and 
one with odd CP. One can form states with orbital angular momenta L = (CP even), 
L = 1 (CP odd), and L = 2 (CP even). 

Alternatively, one can decompose the decay amplitude A into three independent 
components [1^, corresponding to linear polarization states of the vector mesons which 
are either longitudinal (0), or transverse to their directions of motion and parallel (||) 
or perpendicular (±) to one another. The states and || are P-even, while the state _L 
is P-odd. Since J/ip and are both C-odd eigenstates, the properties under P are the 
same as those under CP. 

Consider the polarization three- vectors ej/^ and in the J/ip rest frame. The 
independent decay amplitudes are the rotationally invariant quantities linear in e}^^ and 
60 and involving possible powers of p, a unit vector in the direction of the momentum 
of (f) in the J/ip rest frame. 

The two CP-even decay amplitudes are the combinations e}^^ ■ (contributing to 
Aq and y4||) and e}^^ ■ p ■ p = e*Jj^e*^ (contributing only to Aq), where = p ■ e. 
Equivalently, one can subtract off the longitudinal component of the polarization vectors 
to replace e}^^ ■ by e*Jj^ ■ e'^ , contributing only to A\\, where the superscripts T refer 
to projections perpendicular to p. The CP-odd amplitude e}y^ x ■ p contributes only 
to A±. The case of transverse (|| or ±) polarization states is reminiscent of photon polar- 



ization correlations [ITj] in neutral pion decay. Thus we may write the decay amplitude 
as 



AiB, ^ J/ijct>) = Aoim^/E^)e*j%ef - Aye}}^^ ■ ef /v^ - ^A^e}/^ x e; ■ p/v^ , (6) 

where is the energy of the in the J/ip rest frame, and the individual amplitudes 
are real in the absence of final-state interactions. The amplitudes for the corresponding 
decays of Bg = CPiBg) are Aq = Aq, A\\ = A\\, and A± = —A±. (We can see directly by 
counting powers of p that Aq and A\\ are P-even while A±_ is P-odd.) We have normalized 
the partial widths for the three independent polarization states in such a way that 

dTiB,^ J/ij(P)/dt=\Ao\^ + \Au\^ + \A^\^ , (7) 
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and we may identify 



p{t) = \Ao\^ + \Au 



m(t) = \A 



For purposes of comparing with other notations, we can express the hehcity amphtudes 
Ax (where A = 1, 0, — 1 is the projection of the angular momentum on the x axis) 
in terms of the hnear polarization basis by A±i = {A\\ ± y4_L)/v^, with Aq the same in 
either basis, and in terms of S-, P-, and D-wave amplitudes by 



^±1 = 

With these normalizations 



-D 



An 



-S 



'-D 



dT{Bs J/ij(f))/dt = |Ao|2 + + \A_ 



and 



All 



A, 



\S\^ + \P\^ + \D\ 



P 



V 3 V 3 

The longitudinal and transverse partial widths are given, respectively, by 

dTo/dt = lAol"^ , dTr/dt = \Ai\^ + \A^i\'^ . 
In terms of partial-wave amplitudes, one has 



dt 



1/3S+ y'2/3D|^ 
dVw 



dt 



dt 



2/3^+ Jl/3D| 



2/35 + y'i/3£'r 



|P| 



dt 



while 





2 


m{t) 


P 


2 


p{t)+m{t) |5|2 + |P| 


2 + 


\D 


2 ' p{t)+m{t) \S\^ + \P 


2 + 


D|2 



Finally, we note that in the covariant expression |T 



A 







mim2 



-P2P1 



ic 



17111712 



(9) 
(10) 

(11) 

(12) 



(13) 



(14) 



(15) 



for the decay B V1V2, where e^^^s = and Vi and V2 are vector mesons with masses 
mi and m2 and four- momenta pi and p2, the helicity amplitudes are 



A±i = a ± cVa;2 — 1 , Aq = —ax — b{x^ — 1) 
where x = pi ■ P2/ {'^^1^2)- We thus identify 



(16) 



a{2 + x) +h{x^ - I) 



P = c./2fs2 - V 



D 



a{l - x) - b{x^ - 1) 



;i7) 



Note that S and D both involve a and b. 

The derivation of Eq. is elementary. The is coupled to K^K~ through an 
amphtude ■ {pk+ — Pk-), where the quantities denote 4- vectors. Thus the plane of 
(linear) polarization is related to that of the K^K~ system in the J/ip rest frame. 
By definition, we have taken the linear polarization vector to lie in the x — y plane. 
We may define an angle if) as that of the in the rest frame relative to the helicity 
axis (the negative of the direction of the J/ip in that frame). The spatial components of 
the and J / ip polarizations must be correlated since the decaying strange B is spinless. 
The J/ip then has a single linear polarization state e for each amplitude: In the J/ip 
rest frame, 

A\\ : e = y ; Aq : e = x ; A± : e = z . (18) 
A unit vector n in the direction of the i'^ in J/ ip decay may be defined to have components 

{n^, Uy, riz) = (sin 9 cos ip, sin 9 sin (p, cos 9) (19) 

where (p is the angle between the projection of the on the K^K~ plane in the J/ip 
rest frame and the x axis. The sum over lepton polarizations then leads to a tensor in 
the J/ip rest frame with spatial components (in the limit of zero lepton mass, assumed 
here) 

^ [u-iiv]*[u-ijv\ ~ Lij = 5ij - niUj . (20) 

^±pol 

Physically this tensor simply expresses the fact that massless lepton pairs couple only 
to transverse polarization states of the J/ip, as expected from the structure of the elec- 
tromagnetic interactions. 

Taking account of the definition (^, we then find that the probability for the decay 
Bs {i~^i~)j/^{K~^ K~)^ is proportional to 

^ \A\^ = A,A*U, , (21) 

£±pol 

where 

Ai = AqSix cos 'ip — A\\5iySm'ip / \/2 + iA^5izSi\i^l) / \/2 . (22) 
Consequently, when we use the definitions (|T^, we find 

1 a , , — Tlu = ^ 2|ylo| cos ^ 1 - sm ^cos (^) 

a cos 9 dip acosip at Sztt 

+ sm^ilj{\A\i\'^{l - sin^6'sinV) + \A±\'^ sin^ 9 - Im(A|y4_L) sin26'sinv3} 

+ ^sin2V'{Re(A*A||)sin2^sin2(^ + Im(A*A^)sin2^cos(^}] . (23) 
v2 

The overall normalization has been chosen to agree with our previous conventions when 
one integrates over angles. For Bg decays the interference terms involving A± amplitudes 
are of opposite sign and all other terms are unchanged. 
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Integration over cos ip leads to the distribution 

[\AoWl - sin^^cosV) + IM'^i'^ - sin^^sin^y?) 



dcos9 dip dt Sir 

+ lA^psin^^ - Im (A|Ax)sin2^sin(^] . (24) 

Performing the integrals over if and taking account of the differing time-dependences of 
the decays of -B^^-*, we obtain the result (0). This is a suitable single-angle distribution to 
employ if one wishes to disentangle the CP-even and CP-odd components of the Bg. The 
two- angle distribution ( p^ allows one to separate out the individual quantities \Aq\'^, 
and \A^\^. 

An interesting oscillation appears in the interference terms between CP-even and CP- 
odd decays. For example, in the two-angle distribution (^4]), since the respective time- 
dependences of the B^'^^ and B^^^ decay amplitudes are e"*™-^'*"^^*/^ and ^ 
the term — Im (A|y4_L) behaves as |y4||(0)y4^(0)| sin(Amt — 5)e~'"*, where 5 is a strong 
final-state phase shift difference: v4||(0)*A^(0) = \A\\{<d)A^{<d)\e^\ Thus, if one tags 
the flavor of the decaying Bg, one can observe the effects of Am in the J/ipcf) final 
state through the interference of final states of opposite CP (if both are present). The 
oscillation term averages out to zero if the initial numbers of Bg and Bg are equal. 

The distributions (|), (|23|) , and (p^ also permit one to separate out the components 
|y4oP, l^iiP, and the interference terms Im(y4|y4_L), Re(y4oy4||), and lm{AlA^_) for 

the decays B*^ J/ipK*^. (Here and subsequently we imply the sum over a process and 
its charge-conjugate.) Moreover, in the limit of flavor SU(3) symmetry, one expects the 
ratios of the relative components in B^ J/ipK*^ to be the same as those at proper 
time t = in the decays Bg J /i^fp |[T^- Thus, an analysis of B^ J/ipK*^ can 
provide an independent estimate of the relative contributions of CP-even and CP-odd 
final states at t = to the decays Bg — > J/ipip-, enhancing the ability to determine 
and Pi,. 

If the K* is observed to decay to the CP eigenstate -ft'svr^, the amphtudes Aq and 
^11 refer (as in the case of J/ip4') to the CP-even eigenstate, while A± refers to the CP- 
odd eigenstate. The expected dominance of the CP-even eigenstate (see below) means 
that in B^ J/ipRsir^ events the CP asymmetry will tend to be opposite to that in 
5° J/ipKs ||12|, 0. Since the rates for observing the processes B^ J/ipKs and 
B^ — s> J/ipK*^ — > J/ipRgTi^ are comparable (taking account of branching ratios and 
typical detection efficiencies), the incorporation of J/ipRsir^ data may add statistical 
power to any experiment studying the J/ipKs final state, even when the vr" is not 
observed directly but its existence inferred. 

The distribution (^ permits one to separate amplitudes of opposite parity from one 
another even if the K'^K' system in Bg — >■ J / %l)K^K~ or the Ktt system m B ^ J / ipKn 
is not a vector meson [|l^. This is easily seen by considering the density matrix pij of 
the J/ip, expressed in terms of linear polarization states, so that the decay rate is 
proportional to PijLij, with defined in (|^). If we integrate over (p, we find 

, , /l + cos2^\ . . 

[Pxx + Pyy) 7, + Pzz Sm , (25) 



dcosO dt 
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where pxx and pyy correspond to linear J/%Ij polarization states in the plane of the two 
pseudoscalar mesons, while pzz corresponds to J/ip polarization perpendicular to this 
plane, and thus represents an amplitude with parity opposite to those contributing to 
pxx and Pyy. For cases where each particle in the final state is a C eigenstate as in 

— >• J/ifjK*^ — > J/ipKs'K^ and Bg — >■ J/^Jct) ~^ J/i^KsKi the parity separation is also 
a CP separation and the transversity analysis can be used without the need to extract 
the vector resonance from nonresonant or other background. 

The dominance of the |^oP contribution in B^ J/ipK*^ decays 0, ^ ^ implies 
via flavor SU(3) that the |^oP contribution should also dominate Bg J/'4"P) and hence 
that -B^"-* — > J/ ip'P is likely to be suppressed in comparison with J/i)(j). Thus the 

initial angular distribution is very likely to be dominated by the l + cos^ 9 component. As 
time increases, the fraction of the angular distribution proportional to this component 
will decrease while that proportional to sin^ will increase. It should be possible to 
separate out the two components by a combined analysis in 9 and proper decay time. If 
the sin^ 6 component does not show up even at large times, a single-exponential fit to 
the decay should provide a good estimate of the lifetime of the CP-even eigenstate. 

The angular distribution in (|^) has the form of an ellipsoid which is prolate if p{t) > 
2m{t) and oblate if p{t) < 2m{t). (Here we imagine an average over ip to have been 
performed.) If all three partial waves are equally populated pit) = 2m(t) since there are 
two partial waves with even CP and only one with odd CP. For this case the angular 
distribution is isotropic in cos 9 as expected. 

If p{0) > 2m(0) (i.e., if the CP-even decay is initially more than 2/3 dominant), if the 
CP-even eigenstate -B^^-* has the greater decay rate as expected, and if there is a non- 
zero odd-CP component m(0) ^ 0, then the angular distribution in transversity will be 
initially prolate but will eventually become oblate as the quantity p{t) — 2m{t) changes 
sign. This effect can be noted by dividing the events into two bins with | cos6'| < 1/2 
and I cos 6*1 > 1/2, denoted by E (equatorial) and P (polar) respectively: 

- E = ^ W) - 2m(«)] (28) 

The difference between the numbers in the two bins provides an experimental number 
whose sign will change in time under the assumptions noted above. This two-bin analysis 
of data can be adjusted for optimum statistics by changing the sizes of the bins to 
correspond to the initial values of p{0) and m(0). 

The analysis performed by CDF [|l^ for Bs J/'4"P was also based on a single- 
angle distribution, but it separated out the transverse component from the longitu- 
dinal component. In the absence of vertex cuts these would be, respectively, Tt = 
J^dt{\A\\\^ + l^xH and Fq = /o°° dtdAoP). With a minimum vertex cut of 50 /xm, the 









Jl/2 




P 
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result obtained was ro/(ro + Tt) = 0.56 ± 0.21 (stat) (sys). The transverse com- 
ponent contains both CP-even and CP-odd contributions, while the longitudinal com- 
ponent is CP-even. The separation of transverse and longitudinal components makes 
sense only if the B^'"^ J/'^4' decay is negligible, or if the lifetime difference between 
and can be ignored. 

Corresponding determinations of ro/(ro + Tt) for the decay J/tpK*^ are 

0.65 ± 0.10 ± 0.04 (CDF) 0, 0.97 ± 0.16 ± 0.15 (ARGUS) 0.80 ± 0.08 ± 0.05 
(CLEO) [H, and 0.74 ± 0.07 (world average) |1. This last value is compatible with 
the corresponding one for Bg — J/ip4'- A discrepancy would have indicated either a 
violation of SU(3) or the lifetime effect mentioned above. 

The presence of two eigenstates with possibly differing lifetimes can affect any de- 
termination of t{Bs). When observing the Bg in a final state of definite flavor, such as 
Dg-K or Dgii^e, one will be effectively measuring the average lifetime f of the CP-even 
and CP-odd states. Most measurements reported up to now, including a recent CDF 
determination p3[ leading to a world average of t{Bs) = 1.58 ± 0.10 ps, are of 
this type. This quantity is expected to be very close to the B^ lifetime, for which 
the world average is t{B^) = 1.57 ± 0.05 ps. However, minimum- lifetime cuts can 
bias the sample against the CP-even (shorter-lived) component, leading to results which 
depend on the cut if a single-exponential fit is adopted. 

To summarize, we have found that a combined analysis with respect to proper decay 
time and a single transversity angle in the decay Bs — >■ J /ipcf) can determine the lifetime 
of at least the CP-even and possibly the CP-odd mass eigenstates of the B^ — Bs system. 
Additional information about the properties of the J/ mode at proper time t = can 
be obtained by a similar analysis of the decays B^ J/ipK*^. These analyses can 
already be attempted with the data sample now in hand. 

We thank B. Kayser, R. Kutschke, K. Ohl, M. Paulini, E. A. Paschos, M. P. Schmidt, 
and W. Wester for useful discussions. This work was supported in part by the United 
States Department of Energy under Contracts No. DE AC02 76CH03000, DE FG02 
90ER40560 and W-31-109-ENG-38. 
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